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We introduce chiral-even and chiral-odd meson wave functions as vacuum-to-meson matrix ele- 
ments of bilocal quark operators with well-defined (geometric) twist. Thereby, we achieve a Lorentz 
invariant classification of these distributions which differ from the conventional ones by explicitly 
taking into account the trace terms. The relations between conventional and new wave functions 
are given. 

O . I. INTRODUCTION 

o ■ 

' Recently, we have introduced a group theoretical procedure to decompose local and nonlocal operators [02] , which 

'"JT , are important in different hard scattering processes, into operators of definite twist. This procedure is based on 
[j , • the notion of geometric twist — mass dimension - (Lorentz) spin, r = d — j, originally introduced by Gross and 
Treiman O]. We are interested on nonlocal light-cone (LC) operators and their matrix elements which describe different 
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phenomenological distribution amplitudes, e.g., parton distribution functions as well as hadronic wave functions. The 
classification of these functions suffers from the fact, that one and the same operator by its twist decomposition 
QQ contributes to different wave functions. It is thus necessary to disentangle these various contributions of the original 
^ ' operator and to specify the contributions of the trace terms. 

T— I , Jaffe and Ji Q proposed the notion of dynamical twist (i) by counting powers Q^"* which is directly related to the 

^— I ■ power by which the corresponding distributions contribute to the scattering amplitudes. They classified the parton 

^^ ', distribution functions which correspond to the independent tensor structure of the matrix elements of bilocal quark- 

antiquark operators. Recently, Ball et al. have used this pattern for the classification of leading and higher twist wave 

functions of p- vector mesons in QCD ^^. They found eight independent (two-particle) meson wave functions. One 

^^ I key ingredient in their approach was the use of QCD equation of motion in order to obtain integral representations for 

wave functions that are not dynamically independent. Ball et al. have already pointed out that the geometric twist is 

more convenient to discuss higher twist effects on a reliable basis ||^. On the other hand, the pion wave functions are 

earlier investigated in similar way P,0] . Meson wave functions and form factors have been discussed in the framework 

Oh! of local operator product expansion and in the infinite momentum frame by ||l0|-|l5|] and in the framework of nonlocal 

2:^ operator product expansion [^-Q. 

However, such different notion of dynamical twist is only defined for the matrix elements of operators, it is not 
Lorentz invariant and also not simply related to the contributions of (higher) geometric twist. Quit recently, in order 
to show the relation between the different definitions of twist, we have calculated the forward quark distribution 
functions by the help of our bilocal quark operators with definite twist |l9[ . We proved that the two definitions 
C^ ' of twist do not coincide at higher orders and we gave relation between Jaffe and Ji's distribution function and our 
distributions with geometrical twist. 

The aim of this paper is to present the (two-particle) m,eson wave Junctions which are related to the nonlocal 
LC-operators of different geometric twist. In that framework, it is possible to investigate in an unique manner the 
contributions resulting from the traces of the operators having well-defined twist. 

II. TWIST DECOMPOSITION OF BILINEAR QUARK OPERATORS 

In ||l|,H we explained a procedure of twist decomposition for the bilocal quark operators on the light-cone which are 
relevant for the meson wave functions, too: 

Oa{KiX,K,2x) = u{Kix)^aU{KiX, K2x)d{K2x), (1) 

M[al3]{KlX,K2x) = u{KiX)aa(3U {kiX, K2X)d{K2X) , (2) 
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with the path ordered gauge factor along the straight hne connecting the points kix and K2S:: 



U {kiX , K2x) — P exp < ig / dw x^Afj^{wx) > . 



(3) 



The resulting decomposition for the vector and skew tensor operators is: 

OaiKlX, K2S:) = O^'^^iKlS:, K2X) + 0;7^(Kli, K^i) + Oj7''(Kli, K^X) , (4) 

-|(AJlX, K2i) + Mf-|](Kii, K2i) + Mf-| 



M[„^] (Kii, K2i) = Mf^l (Kii, K2i) + Mf^l] (kiX, K2i) + MM (Kii, K2i) (5) 



with 



0*^^(ki5:, K2i) = X^Ofj,{KiX, K2X) = u{Kix){'-fx)U {kiX , K2X)d{K2x) (6) 

0'^\KiS;,K2i)= f dt[da + ^{\nt)x^DJx^'0^{Kitx,K2tx)\^^- (7) 

O^J^{Kii,K2x)= / di('(5^(2;a)-a;''a„-(l + 21ni)x„a^-(lni)2;„a;''n')0^(Kite,K2te)|^^^ (8) 

0^J'^iKii,K2x)^ia di('(l + lni)a^ + ^(lni)a;^nJ0^(Kite,K2te)|^^- (9) 

^["|](«i5-'«2i) = / di{2ta[^(5^j - (1 - t){2xi^dp]d'^ ~ X[J^^p)}x''M[^,]i>i^tx,K2tx)\^^^ (10) 



'0 





— a; 



[„9^]a;l^a''l)}l'V](«:ite,K2te)|^^^ (11) 



Mf-/](/.ix,/.2i) = / dtij^{x[„5j;x''ln-2a;[4<5j;(x9) -a;[^9^])a'^l}Af[^,](«;ite,«2te)|^^_. 

J 



(12) 



Let us remark that the vector and skew tensor operators of twist r, Oa (nix, K2S:) and M^^i, {kix, K2x) , are obtained 
from the original (undecomposed) operators Oa{Kix,K2x) and M[a0]{Kix , K2x) , Eqs. (^ and (g), by the application 
of the corresponding twist projectors (including the i-integrations), Va and 'Pr^o] , defined by Eqs. (^ - (0) and 
(0) - (|l|), respectively: 

0(;)(KlX,K2£) = (Pi")^0^)(Kl£,K2£) (13) 

M[2,(kix, «2£) = (P[3^''1aV])(ki5;, /«2i) (14) 

with 

(pM x7'("'))^=5"">i")''. (15) 

(p(r)^p(.'))g^5..>M[M^]. (16) 

These projection operators contain the corresponding symmetry of the Young pattern. 

Additionally, the twist-2 vector operator, Eq. (0), is related to the corresponding scalar operator, (||). This leads 
to relations of the corresponding meson wave functions, Eqs. (|l8|) and (|2l|) below. 



III. THE p-MESON MATRIX ELEMENTS OF LC-OPERATORS WITH TWIST 

In this section we define the meson wave functions for the bilinear LC-quark operators with definite twist sandwiched 
between the vacuum and the meson state. As usual, the matrix elements of the meson are related to the meson 



momentum P^ and meson polarization vector e}j , respectively, with P^ ~ mj,, e^^^ • e^'*'-' — —1, P ■ e^"^) — 0, nip 
denoting the meson mass. 

Taking meson matrix elements of Eqs. (g) - ( [l2| ) we see that, observing the correct tensor structure by the use of 
Pp, e}i and x^, we may introduce any parametrization for the matrix elements of the undccomposed operators, e.g., 



(0|O„(^i5;,«2i)|p(P,A)) ^/pm^y du(^el^^^<^iiu,fi^)+P^^^^2{u,fi^)+ic.^ft7^'^3iu,ti^))e'^^<^^P\ (17) 



iP 



[iPf 



where y? denotes the renormalization scale and fp is the vector meson decay constant. We defined k — (ki — K2)/2. 
According to the above projection properties we are able to introduce one (and only one) wave function for any 
operator of definite twist. 

We start with the chiral-even scalar operator. The matrix element of this nonlocal twist-2 operator, Eq. (y), taken 
between the vacuum (0| and the meson state \p[P, A)) reads 

(0|O*-2(«,i,«2i)|p(P,A)) =/pm,(e(^)x) /' du<i>(^)(", A^^) e'^^^^^^ ^ /pmp(e(^)x) f] ^^''^^f ^^" $(f)(^2). (18) 



n=0 



Here, <J>(^^(m, y?) is the twist-2 meson wave function which is related to the corresponding moments 

Jo 
For brevity, we use the shorthand notation 

^ = u - (1 - It) = 2u - 1. 



(19) 



The wave functions describe the probability amplitudes to find the p-meson in a state with the minimal number of 
constituents, e.g., quark which carries momentum fractions u and antiquark with 1 — u, respectively. 

Now we consider the chiral-even vector operator. Using the projection properties (|l3|), together with (|l5|), we 
introduce the meson wave functions ^^'^' {u, fi^) of twist r by 

{Q\Oi^Hnii,n2iMP,X)) ^ (0|(pM^oJ,^^)(«ii, ^2i)|p(P, A)) - /pm,pM^(eJ,^) I dzi$(-)(^,, ^2) gi«5(sp)^ ^ (20) 

which, for t = 2, is consistent with (fsl). Using the twist projection operators as they are determined by Eqs. (R) - 
(0) we obtain for the twist-2 operatorffrom now on we suppress /i^) 



0'J\Kix,K2iMPX))=fpmp dt 



d^ + -{\nt)x^D\{e^^^x) / du $(2)(y) e-«*(-^)|^^^ 



fpuip f dt f du^^^\u)\e'^^ +iK^tPa{e^^'>i) + ^{e^^^S:)ml{iK^t)^{lnt) 
Jo Jo '- 2 



,j.K,e,t{xP) 



n=0 



and for the higher twist operators 
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{Q\0^;^^{kiX,K2x)\p{P,\)) = fpTUp dt dw$(3)(M) ('e(,-^)(£P)-P„(e(^)i)')i<t-Xc,m2(e(^)i)(i<t)2lnt 



^-^ n\ ln+ 1 \ XI 



(^)iN n(n-l)^ m2(e(^)i) 



xPy^(n + 1)2"" {iPf I' 
(0|O;7''(Kii,K2i)|p(P,A)) = i/pmpi«(e(^)i)m2 f dt f du^''^Hu){iKitf{\nt)e"'^'^^P"> 



Jo 



^ [HxPyyi ,,^ n{n - 1) . m^(eWi) 

~Jp^P 2_^ . "i^ 



n=2 



i\ " 2(n+l)2 " (.tP) 
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(21) 
(22) 
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In the first line of any equation we have given the nonlocal matrix element of definite twist, and in the second 
line, after expanding the nonlocal expression, we introduced the moments of the meson wave functions; thereby, the 
^-integrations contribute the additional n-dependent factors. Obviously, the trace terms which have been explicitly 
subtracted are proportional to m^. According to the terminology of Jaffe and Ji, they contribute to dynamical twist-4. 

For the twist-2 operator we observe that the terms proportional to Ca , Pa and Xa have contributions starting with 
the zeroth, first and second moment, respectively. The twist-3 operator starts with the first moment, and the twist-4 
operator starts with the second moment. Analogous statements also hold for the twist-r operators below. 
Putting together the different twist contributions we obtain, after replacing iK^t^xP) by td/dt and performing partial 
integrations, the following matrix element of the original operator (C — k{xP)) 



\OaiKlX, K2x)\p{P, \)) = fpTUp 



e(A)5 /-l 



iP 



du 



($(2)(^)_$(3)(^))[eo(iCC)-ei(iCC)] 



(27) 



+ ei-^) j du(<J>(2)(^)ei(iC0 + '^^^\u)[eo(m " ^MO] 



1 



2""(iP)- .0 
where we used the following "truncated exponentials" 



^cT^ml I du(<i>(2)(^)-2$(3)(y) + $(4)(y)) eo(i«)-3ei(iCO+2 / dteiim) 







eo(ia) = e'«, ei(iCO 



die'«* 



3'« - 1 
iC^ ' 



,e„+i(iCO 



(-1)^ 



dii"e"^«*. 



(28) 



As it should be the application of the projection operators Va onto ( p7| ) reproduces the matrix elements ( |2l| ) - (|2£ 
In comparison with (jlj) we also observe that the wave functions are accompanied not simply by the exponentials, 
eo(iC<^), but by more involved combinations whose series expansion directly leads to the representations with the help 
of moments. Let us mention that the twist-2 part of Eq. (|27| ) and the local expression (|22| ) are in agreement with 
StoU's result |20| for the light-cone twist-2 operator. 

Now we consider the chiral-even vector operator O^aiKix, K2x) — u{Kix)"fa"f5U(Kix, K2x)d{K2x), which obeys re- 
lations Eqs. (^ - (^), as well as (13) and ( p^ with the same projection operator as the vector operator. Let us 
introduce the corresponding meson wave functions S''^-'(z) of twist t by 



|0(;)(«ix,«2i)|p(P,A)) = i(/,-/J^^^^^ duSW(u)e'««^^)) , (29) 



where fj denotes the tensor decay constant. The vacuum-to-mcson matrix elements of these vector operators of twist 



T are obtained as follows: 
{0\OZ\'^ii,^2i)\p{P,X)) = l(fp-fJ^^^^^^^)mpef'^-ef^P,i. [\t / dMS(3)(^)(l -|- i<t(iP))e-«'(^^) (30) 

^ ^ fTT-p ' Jo Jo 



— 2 \jp Jp 



^rriu + md 



I'P'^a 



^13 



Ppx^ 2_^ 



{iK{iP)r ^(3) 



(31) 



n=0 



We see that only the twist-3 operator give a contribution and the twist-2 and twist-4 operator as well as all trace 
terms vanish. Thus, the matrix element of the original operator reads 

{0\O,^{K,i,^2iMP,X)) ^ i(/p - /j "'" + "'" )mp6/^-eyV^x. [' duE^'Hu)eo{ia)- (32) 

z \ Trip / Jo 



The matrix element of the simplest bilocal scalar operator arises as 

{Q\u{K^x)U{K^x,K2i)d{K2i)\p{P,\)) = -i{.fj - fp ^""^ "''' ) (e(^)i)mg / duT^^\u)e'^^^^' 



. I rT r ™f "I" ™'i ^ /' (A) -\ 2 Y^ 

^-nfp -fp^^ — j(« '^)^p2^ 

P n=0 



n! 



(33) 

(34) 



where T'^'(w) is another spin-independent twist-3 meson wave function. 

Now, we consider the matrix elements of the chiral-odd skew tensor operators. The corresponding wave functions 
are introduced by 

{0\Mi^J{n,i,n,xMP,X))^ifJpl:^lf''\{eWp,~e(^^P^ T d^ *M(^) e"^«(s^)) . (35) 

•^ 

The matrix elements of the skew tensor operators of twist r are obtained using the projectors determined by Eqs. (llQ) 
-dll): 

{0\Ml:^^{KiX,K2iMP,X))=ifJ f dt f du^^'Hu)\2te\^^Pp]{2 + iK(tiiP)) 

Jo Jo '- 

+ (1 - t)m2i[„{4(i<i)eJJ' + iiK(t)^{e^^{iP) + P^, (e(^)i)) }] e-«*(*^) (36) 

-'h 2^ „, ^n r^[«^/31+(„ + 2)(n + l)5p''["'''3] +(„ + 2)(n + l)5p'^[n'/31 + ^'^l 5P J/' ^^^^ 
(0|MMl(«iS.«2i)lp(KA))=i/J y 'it^-^f cki<i*m(u)mJi|„{e™+i<tPa,(e»>i)}e"«'<""> (38) 

^ ■'^^ n! 71 + 2 xF ^ \ f^\ '^ xF / 

n—l 

{Q\Ml:^^{K^~x,K2~x)\p{P,X)) = -i/jy''dti^y''d«(ii<)2*(4)(«)™2ij^|e(J)(£P) _P^](eWx)}e-«*(*^) (40) 

i/p 2^ n! " n+liP ["r^] ^'^ iP )■ ^ ' 



Again, the moments of wave functions of twist t — 2,3 and 4 start with n = 0, 1 and 2, respectively. In addition, we 
remark that only those terms of the operator (O) contribute to the twist-3 wave function which result from the trace 
terms of_(M). Analogous to the vector case the forward matrix element of the 'true' twist-3 part of (O) vanishes. 
In Eq. ([40|) only the twist-4 operator contributes which result from the trace terms of (|l^). After multiplication of 
( [40| ) with Xa (or xp) the matrix element vanishes because the corresponding vector operator does not contain any 
twist-four contribution. 

The matrix element of the original skew tensor operator is obtained as 

{0\M[^f,]{K,i,>,2iMP,X)) =ifj\2e\^^Pp] I du^^'\u)eo{m (42) 

[xPf 



■^[c.Pp] (^P)2 j^ d,i{vI'(2)(y)[eo(iCO + 2ei(i«) +6e2(iCO] 



, m? '■^ 



*(3)(u) [1 + 2eo(iCO + 6e2(iCe)] + ^^'H") [l + ^oiKO - 2ei(i«)] } 
+ ^["4]^^y du{vI'(2)(^,)[eo(iCe) -2ei(ia) -2e2(iC0] 

*(3)(u)[i + 2e2(iCO] - *''Hw)[l + eo(i«) - 2ei(iC0] }" • 



Now we finish the determination of the eight meson wave functions which result from the nonlocal light -cone quark 
operators of definite twist. 

IV. RELATIONS BETWEEN NEW AND CONVENTIONAL MESON WAVE FUNCTIONS 

Obviously, since these new wave functions of mesons are related to true traceless operators they differ from the 
conventional ones |^,0| for higher twist at least by the contributions from the trace terms. As far as the scalar LC- 
operators are concerned which definitely are of twist-2 the new and the old wave functions coincide. However, for the 
vector and (skew) tensor operators also the contributions of dynamical twist-2 differ from those of geometric twist-2. 



Now, we are able to compare our new and the conventional meson wave functions. We rewrite the matrix elements 
@j, (11) and @ by choosing 



P, 



^Xj fy tt' ^ 



a -^a 



-^((.P)-/^ 



1. mi 



-S^'>x 1 , n e^^^x 



'J' ^""^"+2^"iP' """ """^"^^"^^^ 2'^""'^'(£P)' 



(43) 



.(^) 



where pa is a light-like vector (p^ = and p ■ x = P ■ x) and e^^^j^ is the transversal polarization vector of the p-meson. 
Our result is 



\Oa{KlX, K2x)\p{P, X)) = fpTTlp 

1 , m^{e'-^'>x) '■1 






(2)(i,)eo(i«) +eV'] / d^.{<i>(2)(^.)el(i«) + $^'H«)[eo(i«) -ei(iCC)]} 



2"" (iP)2 



du 



{<i>(4)(M)[eo(iCO -3ei(iC0 +2 / diei(iat)] - <i>(')(u) [ei(iCO - 2 / dtei(iC^i) 



+ 4$(3)[ei(iCa-y"diei(iat)]} 
(0|O5J«:i5:,/.25)|p(P,A)) = l(fp- fj!!h^±^)mpeJ^-^e'^llp,S:, f duE'^'Hu)eo{m, 



(0|M[„^](^ix,K2i)|p(P,A))-i/J 2eV7„P0] / duvI;(2)(y)eo(i^^) 

+ 2x[aP/3] (^p)2 y^ d7.|2vI/(2)(„)e2(iCe) - vI/(3)(«) [eo(iCe) + 262(1^)] } 

2 ,.1 



(44) 
(45) 
(46) 



*["4uip / d«{2*('^(^)h(iCe)+e2(iCe)] -*^'^(^)[l + 2e2(iC6] +*^'^(")[l + eo(iCe)-2ei(iCO]} 



Comparing these expressions with the meson wave functions with dynamical twist given by Ball et al. |6|J3] we observe 
that it is necessary to re-express the truncated exponentials and perform appropriate variable transformations. After 
such manipulations we obtain the following relations, which allow us to reveal the interrelations between the different 
twist definitions of meson wave functions: (For simplicity, we give here only the integral relations of the meson wave 
functions in a state with quarks which carry momentum fraction u.) 



011 (w) 


= $(2)(y) 


9± (u) 


= <i>(3)(„) 


93{u) 


= $(4)(^) 




= S(3)(u) 


h^iu) 


= T(3)(«) 


4>Au) 


= Vl/(2)(„) 


hf{u) 


= *(3)(^) 


hsiu) 


^¥^\u) 



' ^|($(2) _ 4ci>(3) + 3<i>(4)) (^) + 2 In (-) ($(2) _ 2$(3) + $(4)) („) 1^ 



2u /' J(vI/(2)_vI/(3))(„)^ 



(47) 
(48) 

(49) 

(50) 
(51) 
(52) 

(53) 



' ^|2(*(2) _ M,(4)) («) _ 2!^(vI/(2) _ vl;(3)) (^) _ s[-) (vI/(3) _ vl/(4)) („) |. (54) 

We observe that both decompositions coincide in the leading terms, but differ at higher twist. For instance, the meson 
wave functions g^ (u) and hi (u) with dynamical twist t = 3 contain contributions with geometrical twist t — 2 and 
3. Additionally, dynamical twist t — A meson wave functions 53 (u) and h^^u) contain contributions with geometrical 
twist T = 2, 3 as well as 4. 

Additionally, the conventional wave functions can be written in terms of our new functions. The nontrivial relations 



<I>(4)(i.)-.g3(i.) + 4 / di;^353-4gi^^+0||)(«) + ^ / dvv(l ~ -){g, ~ Ac/l^ +3^\\) (v) , (56) 

2 /•! 



'J'(3)(u) = /i[*)(u) + - / dv{h\''> - 0^) (v) , (57) 

*(*H«) = /i3(ti) + 4 / d«f(/i3-/j[|')(t^)-4 / dt-i;(l--)(/i[*^-0i)(t;). (58) 

The relation between the moments may be read off from Eqs. ( fji] ) - (^) as foUows: 

0||n-*l'\ (59) 

53„ = $i*' - ^ ($i^) - 4<i>(f ) + 3<i>(4)) + j-^ ($1^) - 2$(3) + $(4)) , (61) 

ffS = 2i3\ (62) 

/^£ = T(f), (63) 

0±n = *i'\ (64) 

/^3„ = *1^) + ^ (*1^) - *(^)) - -^ {^l? - ¥,?) 5„o (*i^) - ^i-^)) . (66) 

In terms of the moments the relations between old and new wave functions may be easily inverted; for the wave 
functions itself the expression of the new wave functions through the old ones is more involved. The inverse relations 
are: 

'^'^'-9^ll + ll{9^ll-hn), n>Q (67) 

<^i'^ = 93n + ^ (3.g3„ - 45^:',) + 0[|„) + ^^^\ ^^ (g3„ - ^g^ll + 30||„) , n > 1 (68) 

M.(4) = /,3„ + _1_ (/,3„ _ /,(;)) „ _^ (;,(*^ _ 0^„) , „ > 1. (70) 

V. CONCLUSIONS 

In this letter, we have discussed the calculation of the vacuum-to-meson matrix elements for nonlocal LC-operators 
using the notion of geometric twist. We have found eight meson wave functions with geometric twist t. From the 
field theoretical point of view this Lorentz invariant classification is the most appropriate frame of introducing wave 
functions since the separation of different (geometric) twist is unique and independent from the special kinematics 
of the process. An important result of our calculations are the relations between the new meson wave functions and 
those given by Ball et al. KM- These integral relations reveal the connection between the geometric and dynamical 
twist definitions for the meson wave functions. 

An advantage in our approach was that we use operators with well-defined twist. Therefore, we have not used 
operator relations for different geometric twist in order to isolate contributions of (geometric) twist-3 and 4 (without 
gluonic contributions and neglecting quark masses) what is used in PJ20[|. 
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